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ABSTRACT 

We have obtained high-resolution spectroscopy of ten red giants in the Carina dwarf spheroidal (dSph) 
galaxy with the UVES spectrograph at the ESO Very Large Telescope in order to study the detailed 
chemical evolution of this Galactic satellite. Here we present the abundances of O, Na, Mg, Si, Ca, 
Ti and Fe. By comparison of the derived iron abundances [Fe/H] with metallicities based on the well 
established calcium triplet (CaT) calibration, [Fe/H]caT, we show that the empirical CaT technique 
yields good agreement with the high-resolution data for [Fe/H] > —2 dex, but tends to deviate from 
these data at lower metallicities. With [Fe/H]^ —1.7 dex the mean iron abundance of our targets is 
fully consistent with the peak metallicity of Carina as derived from medium-resolution spectroscopy and 
previous photometric studies, all calibrated onto iron via Galactic globular cluster scales. We identify 
two metal poor stars with iron abundances of —2.72 and —2.50 dex. These stars are found to have 
enhanced [a/Fe] ratios similar to the elemental ratios of stars in the Milky Way halo. In this context 
it is conceivable that the moderately metal poor halo stars may originate from an early dSph accretion 
event. The bulk of the Carina red giants exhibit a depletion in the [a/Fe] abundance ratios with respect 
to the Galactic halo at a given metallicity. One of our targets with a moderately low [Fe/H] of —1.5 dex 
is considerably depleted in almost all of the a-elements by ~ 0.5 dex compared to the solar values. Such 
low values of the ratio of a-elements to iron can be produced by stochastical fluctuations in terms of an 
incomplete mixing of single Type la and Type II supernova (SN) events into the interstellar medium. 
Moreover, the system's slow star formation (SF) rate grants sufficient time for SNel to occur. Our 
derived chemical element ratios are consistent with the episodic and extended SF in Carina previously 
derived from analyses of its color-magnitude diagram. We find a considerable star-to-star scatter in the 
abundance ratios. This suggests that Carina's SF history varies with position within the galaxy, with 
incomplete mixing. In addition, or alternatively, the SF rate is so low that the high-mass stellar IMF is 
sparsely populated, as expected on statistical grounds in low-mass star clusters, leading to true scatter 
in the resultant mass-integrated yields. Both ideas are consistent with slow stochastic SF in dissolving 
associations or star clusters, so that one may not speak prima facie of a single "SF history" at a detailed 
level. 

Subject headings: Galaxies: abundances — Galaxies: dwarf — Galaxies: evolution — Galaxies: stellar 
content — Galaxies: structure — Galaxies: individual (Carina) — Local Group 



1. INTRODUCTION 

Dwarf spheroidal galaxies (dSphs) are intriguing ob- 
jects: they are the least luminous, least massive galaxies 
known (e.g., Grebel et al. 2003). They also appear to be 
the most dark-matter dominated objects (e.g., Gilmore et 
al. 2007 and references therein). For many years dSphs 
were considered the potential building blocks of larger 
galaxies, and thus potential key stones from which to ex- 
tract the merging history of the Milky Way (e.g., Gallagher 
& Wyse 1994 and references therein). While there is clear 
evidence for past mergers of dwarf galaxies with the Milky 
Way, the age distributions and metallicity distributions of 



stars in the dSphs, plus recent detailed abundance anal- 
yses (Grebel & Gallagher 2004; Venn et al 2004), have 
called a simple building block scenario in question. These 
limit accretion and merging of systems like the surviving 
dSphs to early epochs only (Unavane et al. 1996). 

Among the dSphs orbiting the Milky Way the Carina 
dSph stands out because of its unusual star formation (SF) 
history. It is the only dSph to exhibit clearly episodic SF 
interrupted by long quiescent periods. Carina has been the 
subject of numerous, primarily photometric, studies (for 
an overview of past work, see Koch et al. 2006a; hereafter 
Paper I) . We recently analysed Carina's metallicity distri- 
bution function (MDF) derived from medium-resolution 
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spectroscopy around the ncar-infrarcd Calcium II triplet 
(CaT) of 437 red giants (Paper I) , from which we found a 
mean metallicity of [Fe/H]caT = —1.7 dex on the scale of 
Carretta fc Gratton (1997), with a wide range from about 
—3.0 to 0.0 dcx. Hence, there is evidence of considerable 
enrichment during the evolutionary life time of this dSph 
from about 12 Gyr ago to a few Gyr ago. 

Low- to medium resolution spectroscopic studies are 
well suited to investigate the true, spectroscopic mean 
metallicity and metallicity spread of individual dSphs and 
to derive their MDF. Combined with photometry, this per- 
mits one to assess to what extent age and metallicity con- 
spire (i.e., the age- metallicity degeneracy in color-space), 
to search for possible spatial gradients, and to explore the 
galactic evolutionary histories taking into account their 
chemical evolution. One should note that the overall 
'metallicity', which is commonly represented via [Fe/H] 
or [Fe/H] CaT, is generally most sensitive to the integrated 
past SF history (see, e.g., Paper I), while it is rather in- 
sensitive to the details of that SF history. Unravelling SF 
histories requires knowledge of individual elemental abun- 
dances and of the associated stellar ages: the evolution of 
specific chemical element ratios is intimately linked to the 
SF history of a galaxy, the initial mass function (IMF), 
and the binary fraction governing the rate of Supernovae 
(SNe) of Type la. 

An important example of such abundance ratios is 
[a /Fe] , which can be efficiently used as a "chemical clock" 
(Tinsley 1976; Wyse & Gilmore 1988; Matteucci 2003). 
The a-elements (O, Mg, Si, Ca, Ti) form in quickly evolv- 
ing high-mass stars, which end their lives in core-collapse 
SNe II. Fe is also synthesized in SNe II, but predominantly 
in SN la, which occur on much longer time scales. Hence, 
modulo loss of metals from the galaxy, the [a/Fej-ratio of 
a a stellar population records its SF history, reflecting the 
relative delay in the injection of Fe and a-elements into the 
interstellar medium (Matteucci & Grcggio 1986; Gilmore 
& Wyse 1998; Lanfranchi & Matteucci 2004; Lanfranchi 
et al. 2006; Robertson et al. 2005; Font et al. 2006). 

Stars in dSphs have lower [a/Fe] abundance ratios on 
average than do Galactic halo stars of the same metallicity 
(Shetrone, Cote & Sargent 2001; Fulbright 2002, Shetrone 
et al. 2003 [hereafter S03]; Tolstoy et al. 2003 [hereafter 
T03]; Geisler et al. 2005). This is generally interpreted 
as evidence of a low SF rate in dSphs compared to the 
Galactic halo, as predicted by Unavane et al. (1996). The 
actual values will sensitively depend on the SF and ISM 
mixing histories, and the IMFs, of both the dSph galaxy 
and of the formation location of the field halo stars. 

From their analysis of five red giants in Carina, S03 
and T03 deduce the imprints of SF bursts and pauses in 
Carina's abundance patterns. Each SF burst can be ex- 
pected to be accompanied by a rapid increase in [a/Fc] 
(Gilmore & Wyse 1991, 1998), followed by a smoother 
decline in this ratio during any subsequent period of low 
SF. This simply reflects the continuing production of Fe 
in SNe la, which continue for many Gyr after the forma- 
tion of their progenitor stars (Matteucci & Greggio 1986), 
with little a-enrichment due to the paucity of SNe II. This 
SF and consequent enrichment history was claimed to be 
also reflected in variations of the heavy-element patterns, 
such as some [s/r]-process element abundance ratios vs. 



metallicity. The illustrative scenarios proposed by S03 
and T03 suggest an ancient epoch of SF (around 13 Gyr 
ago), a subsequent quiescent phase of some 3 Gyr, and 
another SF period approximately 7-11 Gyr ago. This is, 
however, not fully compatible with Carina's detailed SF 
history as derived from color-magnitude diagram (CMD) 
modeling (e.g., Smecker-Hane et al. 1994; Mighell et al. 
1997; Hurley-Keller et al. 1998; Monelh et al. 2003; Rizzi 
et al. 2003) and thus merits further investigation. 

Our current study is a further step into this direc- 
tion. We have enlarged the previously available sample 
by adding ten new stars, for which we obtained high- 
resolution spectra using the UVES spectrograph at the 
Very Large Telescope (VLT) of the European Southern 
Observatory. Considering the considerable expense in ex- 
posure time for high-resolution studies, this represents al- 
ready a substantial improvement, even though these data 
still only yield snapshots in sampling the SF episodes of 
Carina. Each of our target stars provides a valuable pic- 
ture of its local environment at the time of its formation. 
We have underway a much larger VLT/UVES study of 
Carina, to extend this work. 

Calibration of the reduced width of the near-infrared 
CaT index in terms of metallicity is another issue in- 
volving an understanding of the abundance of the a- 
elements. The widely used, empirical calibrations (Ar- 
mandroff & Da Costa 1991; Rutledge et al. 1997a,b) are all 
based on the metallicity scales of Galactic globidar clus- 
ters (e.g., Zinn & West 1984; Carretta & Gratton 1997; 
Kraft & Ivans 2003), themselves calibrated by iron el- 
emental abundances derived from high-resolution echelle 
spectra. Cole et al. (2004) have extended this calibration 
to ages as young as 2.5 Gyr and metallicities in the range 
of —2.0 < [Fe/H] < —0.2 dex so that even the metallicities 
of stars in dSphs with extended SF histories and promi- 
nent intermediate-age populations can be measured from 
the CaT. However, [a/Fe] varies systematically across this 
range in the calibrators, so that it remains an open ques- 
tion whether and to what extent the extrapolation of the 
calibrations towards the metal-poor and metal-rich tails of 
the MDFs introduces a possible bias or significant devia- 
tions from the "true" Fe-line stellar metallicity. The CaT 
calibration explicitly assumes a [Ca/Fe] ratio that is simi- 
lar to the one found in the Galactic globular clusters, i.e., 
a ratio of the order ~ -1-0.3 dex (Carney 1996; Mc William 
1997; Pritzl et al. 2005). However, in dSphs this ratio is 
not an a priori known quantity, and generally this value 
shows a wide range, even within one galaxy. There is ev- 
idence that some stars in some dSphs can reach [Ca/Fe] 
values as enhanced as the Galactic halo (Shetrone et al. 
2001, SOS). Hosier et al. (2007) argue that the CaT line 
strength is naturally best mapped onto [Ca/H], where the 
respective calibration is independent of the SF history of 
the system and valid over a broad range of metallicities 
and ages. Our measurements of [a/Fe] in Carina will then 
aid in elucidating the validity and accuracy of traditional 
CaT metallicity measurements in the dSphs. 

In this current study we focus on the abundance pat- 
terns for Na, Fe and the a-elements O, Mg, Si, Ca, Ti in 
ten red giants in Carina. An analysis of the heavy element 
abundance ratios will be presented in a forthcoming paper 
(Koch et al., in preparation). In §2, we describe our data 
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and their reduction. §3 describes the atomic data and at- 
mospheric parameters used in the abundance analysis in- 
cluding a detailed discussion of errors. In §4 we present our 
abundance results for iron and the individual a-elements 
before discussing their implications for Carina's chemical 
evolution in §5. §6 summarizes our findings. 

2. DATA AND REDUCTION 

Our data were obtained in the course of the ESQ Large 
Programme 171.B-0520(A) (PI: G. F. Gilmore; see also 
Paper I; Wyse et al. 2006; Wilkinson et al. 2006; Koch ct 
al. 2007, for details) that aims at studying the kinematic 
and chemical characteristics of Galactic dSphs. Two of 
these galaxies were extensively studied in medium resolu- 
tion mode, namely Leo II (Koch et al. 2007) and Carina 
(Paper I). In parallel, a number of high-resolution spectra 
of red giants in Carina were obtained, in order to study 
the detailed chemical composition of this particular galaxy. 
To allow us to perform accurate chemical abundance stud- 
ies at a high signal-to- noise (S/N) ratio within a reason- 
able integration time, 30 bright giants below the tip of 
the RGB were selected from the CMD, reaching V-band 
magnitudes slightly above 18 mag (see Table 1). Eleven 
of these stars were also observed with medium-resolution 
around the CaT region, from which we derived Carina's 
MDF in Paper I. For these stars we can compare directly 
the stellar metallicities from the CaT calibration with the 
iron abundance from the high-resolution measurements, 
and thus assess the validity of the CaT calibration method. 

2.1. Data acquisition 

In parallel with the observing runs in December 
2003 (see Paper I) , during which medium-resolution spec- 
troscopy of red giants in Carina was obtained, we used 
the high-resolution spectrograph UVES (Ultraviolet and 
Visual EchcUe Spectrograph) at the ESO/VLT in multi- 
object mode. Using UVES in combination with the 
FLAMES fibre facility permits one to feed nine fibres to 
the UVES instrument, where we dedicated one fibre per 
setup to observing blank sky. This configuration facilitates 
the sky subtraction later on. 

The FLAMES/UVES fibre mode has been designed for 
the red arm of the spectrograph. Our observations used 
the standard setup with a central wavelength of 580 nm, 
which gives a total spectral range of 470-680 nm. With 
a fibre aperture on the sky of 1" and our chosen CCD 
binning of 1 x 1 pixels the instrument provides a resolving 
power of --47000. 

2.2. Data reduction 

The data were reduced using the standard UVES 
pipeline (Modigliani 2004^. which operates within the 
MIDAS environmenlPl. The respective standard steps 
comprise bias subtraction and division by a master slit flat 
field. The definition of the spectral orders is obtained by 
the usual means of processing odd- and even-numbered 

This document is available online via 

|http : / /www. eso . org/pro j ect s/df s/ df s - shared/web/vlt/vlt- instrument-pipelines ■ html | 

The European Southern Observatory Munich Image Data Analysis System (ESO-MIDAS) is developed and maintained by the European 
Southern Observatory. 

^■^ IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science Foundation. 



fibre flat field frames, combined with a first guess or- 
der position table from the Th-Ar wavelength calibration 
lamp. In the next step, the wavelength solution is de- 
termined via calibration exposures from this Th-Ar lamp, 
followed by the extraction of the target (and sky) spec- 
tra using an optimal extraction algorithm. The final ex- 
tracted, wavelength-calibrated spectra with contiguously 
merged orders are then sky-subtracted using the dedicated 
sky spectrum taken in each setup. Finally, the individual 
exposures are shifted to the heliocentric standard of rest 
and coadded, weighted by the spectrum's individual S/N. 
Typical S /N ratios of the final spectra achieved in this way 
are approximately 20pixel~^. 

2.3. Membership estimates 

The likely Carina membership of the target stars was 
assessed by means of their radial velocities. We gener- 
ated a template spectrum from a line list comprising 101 
strong absorption features with equivalent widths (EWs) 
as expected for typical red giant spectra in dSphs. Cross- 
correlation of the observations against this template using 
iraf'^H /a^cor task finally yielded accurate radial velocities 
with a median uncertainty of 0.9kms~^. 

For the 11 targets for which we also obtained medium- 
resolution spectra of the CaT region, the velocity measure- 
ments from the UVES spectra are in excellent agreement 
with the velocities derived from the CaT (see Paper I) to 
within their uncertainties, and we find a median offset be- 
tween both estimates of 0.5 km s^^ with a Icr-scatter of 
l.lkms"^ 

Out of our 30 targets, 12 stars clearly peak at Carina's 
systemic velocity around a mean of 225kms~^ with a dis- 
persion of 6.4kms~^. This is in good agreement with Ca- 
rina's radial velocity distribution (Mateo 1998; Majewski 
et al. 2005; Paper I). The remaining 18 targets have ra- 
dial velocities that deviate by more than 120 km s'^ from 
this mean, which corresponds to 20 a and unambiguously 
rules out their membership in Carina. These stars are not 
considered in the current study and are likely Galactic in- 
terlopers (Wyse et al. 2006). We note that the candidates 
for observation were selected deliberately to extend well 
to the blue and to the red of the Carina RGB, so that 
any very metal-poor or metal-rich Carina stars would be 
included. No such stars were found. Two of the radial 
velocity members were identified as carbon stars via their 
strong molecular bands. These spectra were omitted from 
the abundance analyses presented here. 

3. abundance analysis 

For the abundance determinations, we employed the 
common technique of an EW analysis. For this purpose, 
we used the 2002 version of the stellar abundance code 
MOOG (Sneden 1973), in particular using its abfind driver . 
For the van der Waals line damping, the standard Unsold 
approximation was adopted within MOOG. 

The accurate derivation of stellar abundance ratios sen- 
sitively relies on the choice of proper atomic data, in par- 
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ticular the oscillator strengths. Concern about using solar 
^/-values for the analysis of our late-type red giant tar- 
get stars prompted us to perform a differential abundance 
analysis relative to Arcturus {a Boo), a bright K1.5 disk gi- 
ant of only moderate metal-deficiency ([Fe/H]=— 0.50 dex) 
and with atmospheric parameters very similar to those ex- 
pected for the Carina red giants of our sample (e.g., Peter- 
son et al. 1993; Friel et al. 2003; Fulbright et al. 2006, 2007; 
Koch & Mc William 2008). For this purpose, we used the 
Arcturus spectral atlas of Hinkle et al. (2000) to obtain a 
set of reference measurements. 

3.1. Equivalent widths and gf -values 

In order to identify suitable transitions for our abun- 
dance determinations, we compiled a line list from various 
literature sources. Essentially, we exploited the compila- 
tions from the bulge giant and thick disk star studies of 
McWiUiam et al. (1995) and Prochaska et al. (2000), from 
which we assembled a number of absorption features for a 
number of chemical elements (see references in these works 
for the sources of the respective atomic data). We chose 
the lines with the least amount of stellar and telluric con- 
tamination, which thus appear reasonably unblended. 

The EWs of these lines were measured using IRAF's 
splot-task assuming a single Gaussian line profile. Typical 
splot-eiTOTS on the widths are of the order of 5% (a few 
mA), and the EWs of those lines in common with Ful- 
bright et al. (2006, 2007) agree to within 3% with their 
measurements . 

The oscillator strengths used in the majority of abun- 
dance studies are nowadays based on laboratory measure- 
ments, such as the O'Brian values (O'Brian et al. 1991), 
Hannover measurements (Bard & Kock 1994) or the Ox- 
ford 5/- values (Blackwell et al. 1995), all of which are 
employed in the work of McWiUiam et al. (1995) and 
Prochaska et al. (2000). StiU, as Prochaska et al. (2000) 
note, even with such maximally accurate laboratory mea- 
surements at hand, the ^/-values impose a major error 
source in red giant abundance studies, particularly when 
referencing one's abundances to the solar meteoritic (An- 
ders & Grevesse 1989) or the updated empirical (Grevesse 
& Sauval 1999; Asplund et al. 2005) solar abundance ref- 
erence frame. Therefore we have chosen to conduct our 
analysis differentially with respect to the well-studied red 
giant Arcturus, thus minimizing the uncertainties associ- 
ated with solar and laboratory ^/-values. Since Arcturus 
exhibits atmospheric properties similar to the stars of our 
own study, such an analysis will also aid the reduction of 
systematic errors, such as a priori unknown deficiencies in 
the model atmospheres or the influence of weak, unknown 
blends, which might scale with metallicity (e.g., Fulbright 
et al. 2006). To this end we measured the EWs of the lines 
in Arcturus from the high-quality spectral atlas of Hinkle 
et al. (2000). 

Adopting the well established atmospheric parame- 
ters of aBoo (Toff = 4290K, log5=1.55, ^ = 1.67kms-\ 
[Fe/H] = —0.5) from Fulbright et al. (2006) and our mea- 
sured EWs, we derived the (//-values on a line- by-line ba- 
sis in order to reproduce each element's abundance from 

These atmosphere grids can be downloaded from pittp : //cf aku5 . 
See http : //wwwuser . oat . ts . astro . it/castelli 

The reduced and cahbrated photometry is available from the EIS 



Fulbright et al. (2006, 2007). As a result, our Arcturus- 
based oscillator strengths agree with the laboratory values 
to within ~8% on average. The full line list is given in 
Table 2, together with the respective EW measurements 
for each star. 

3.2. Model Atmospheres and stellar parameters 

Throughout our analyses we interpolated the model at- 
mospheres from the updated grid of R. L. Kurucz's[3 one- 
dimensional 72-layer plane-parallel line-blanketed models 
with the convective overshoot option switched off and as- 
suming that local thermodynamic equilibrium (LTE) holds 
for all species. Moreover, our model calculations incor- 
porated the new opacity distribution functions, ODFNEW, 
provided by F. CastellH (Gastelh & Kurucz 2003) for 
the Carina stars, assuming scaled-solar elemental abun- 
dances, and the respective a-element enhanced models, 
AODFNEW, for Arcturus. The latter are necessary since it 
is known that Arcturus exhibits a considerable enhance- 
ment in these elements (Peterson et al. 1993; Fulbright et 
al. 2007) of [a/Fe]- -1-0.4 dex. 

The initial estimates for the stellar effective tempera- 
tures were derived using dereddened colors from the pho- 
tometry obtained by the ESO Imaging Survey (EI^3; 
Nonino et al. 1999; Paper I), complemented by infrared 
colors, which were drawn from the Point Source Catalogue 
of 2MASS (Cutri 2003). For Carina, we use a redden- 
ing of E(B-V) = 0.06 (Mighell 1997), which is also con- 
sistent with the value we found from the extinction maps 
of Schlegel et al. (1998). To correct the infrared colors, the 
respective reddening relations from Cardelli et al. (1989) 
were applied. Finally, photometric effective temperatures 
were calibrated using the empirical relations derived by 
Alonso et al. (1999, 2001). Since the Alonso et al. cali- 
brations of B— V and V— K require colors in the Johnson- 
Cousins (JC) system, we employed the photometric trans- 
formations described in Paper I in order to convert the in- 
put B and V magnitudes from the photometric system of 
the EIS to the standard JC system. The 2MASS K-band 
magnitude was transformed to the Bessel & Brett (1988) 
homogenized system, using the equations provided by the 
Explanatory Supplement to the 2MASS AU Sky Data Re- 
lease (Cutri 2003). Since this transformation includes two 
steps of calibrations, i.e., from the 2MASS to the CIT sys- 
tem and from CIT to the Bessel & Brett (1998) system, 
the transformation errors from both calibrations have to 
be incorporated to yield the final photometric uncertainty 
of the colors (see also Sivarani et al. 2004). Table 3 lists 
the photometric properties of all target stars with radial 
velocities that are consistent with membership in Carina. 

As initial guesses for the metallicity to be used in the 
Toff calibrations and in the model atmospheres we adopted 
the value derived from our CaT study of Paper I for the 
six member stars in common. Otherwise Carina's mean 
metallicity of —1.7 dex (Paper I) was assumed for the re- 
maining four red giants for which no metallicity measure- 
ment was available. 

In practice, we adopted the Toff (V— K) value as the fi- 
nal stellar photometric temperature, since the large wave- 

cf a . harvard ■ edu/ grids ■ html | 
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length baseline of this color index makes this temperature 
the least susceptible to the underlying photometric uncer- 
tainties, yielding the most accurate results. The temper- 
atures obtained from the B— V colors are systematically 
cooler than their respective V— K counterparts by ^70 K 
on average. The random errors on the effective tempera- 
ture, based on the full photometric uncertainties, amount 
to 70 K on average. 

The stellar surface gravities, log 5, were initially deter- 
mined using the basic stellar structure equations 

a M T ft 

log A = log + 4 log - 0.4 (Mboi, - My - BC), 

9q Mq J-eff,© 

(1) 

where we adopted Tofr,© = 5777K, Mboi, © = 4.74 and 
loggi© = 4.44 and a distance of Carina of 94kpc (Mighell 
1997). Moreover, the photometric temperature as derived 
above was used in this relation, and the stellar mass was 
taken to be 0.8 M©, in correspondence to the typical main- 
sequence turn-off mass. The bolomctric correction (EC) 
was obtained by interpolation of the respective Kurucz- 
atmosphere grids. 

Adopting this set of values as the input parameters, 
the final spectroscopic parameters were determined from 
a subset of the clean Fei lines listed in Table 2, for which 
—5.4 < log;^o(E'W/A) < —4.5, corresponding to a range 
of ~20 to 160 mA at 5000A. The lower limit of this selec- 
tion prevents weak lines being spuriously detected, thus 
artificially raising the mean iron abundance near the cut- 
off. The strongest, saturated lines were pruned from the 
list, since they do not lie on the linear part of the curve of 
growth and are predominantly formed in the highest lay- 
ers, where the atmosphere models become progressively 
unreliable (McWiUiam et al. 1995; Johnson 2002). 

Atmospheric parameters were determined by iteratively 
modifying the models to satisfy three constraints in suc- 
cessive steps. As a first step, we determined the micro- 
turbulent velocity ^. For this purpose abundances from 
the neutral iron lines were iteratively calculated, demand- 
ing that there is no trend of abundance with the reduced 
width EW/A. The scatter around the zero slope from a 
least-squares fit was then taken as a measure of the uncer- 
tainty on the microturbulent parameter. In this way, we 
conservatively obtain typical accuracies of 0.1 kms~^ on ^. 

In the second step, spectroscopic temperatures were ob- 
tained by requiring excitation equilibrium. Based on the 
range of Tcff that produces acceptable fits in the Fe i- 
abundance versus excitation potential (EP) plot, we es- 
timate the typical errors of the spectroscopic temperature 
to be ±100 K, slightly larger than the photometry-based 
errors. Generally, our spectroscopic Toff- values are sys- 
tematically higher than the photometric estimates, and 
we note an offset of -1-250 K on average for our standard 
atmospheres. For some of the targets this discrepancy 
can be as high as 500 K. Although a systematic offset of 
^250 K may raise some concern, differences between the 
spectroscopic temperature scales do not come as a surprise 
(see also the discussion in Johnson 2002). In our analy- 
ses we use the spectroscopic temperatures to define the 
atmospheric Toff for several reasons. First, the photomet- 
ric errors of the infrared colors are rather high. The me- 
dian photometric uncertainty of our targets from 2MASS 
amounts to 0.1 mag in V— K. A reddening uncertainty of 



0.1 mag leads to a Tefj variation of ±100K. It turns out 
that an unreasonably large systematic decrease of 0.4 mag 
on average in V— K (up to 0.71 mag for the most deviant 
stars), be it due to photometric or reddening uncertainties, 
would be required to reconcile the spectroscopic and pho- 
tometric temperatures. However, this would correspond to 
an unreasonably large increase of Carina's distance mod- 
ulus of 0.45-0.8 mag (reflecting a distance shift of ~18- 
30kpc). Inconsistent photometric and spectroscopic Tcff- 
values could occur due to errors in the (^/-values, which 
may correlate with the EP (Johnson 2002). However, this 
has been shown not to be the case either for laboratory os- 
cillator strengths (e.g., Blackwell et al. 1995), nor should 
it be a problem in our differential line-by-line analysis. Fi- 
nally, Johnson et al. (2006) have shown that the use of 
a photometric Tefi yielded highly diverging non-standard 
abundance ratios for their sample of old red giants in LMC 
clusters, again strengthening our confidence in the use of 
the spectroscopic values. 

In the third step, the appropriate surface gravity for 
each star was constrained by assuming ionization equi- 
librium, thus requiring that the log g-scnsitive Fell abun- 
dance matches the abundance from the Fei lines, £(Fei). 
These abundances were repeatedly computed for both 
species with loggi varying in small steps until the results 
from both ionization stages agreed reasonably well. By 
demanding that both ionisation levels do not deviate by 
more than 0.1 dex we could effectively constrain our grav- 
ity estimate to within an accuracy of ±0.2 dex. 

As for the case of Teff, the gravities determined from 
ionization equilibrium are also systematically higher than 
the photometrically derived values (eq. 1). This discrep- 
ancy cannot be resolved by adopting our spectroscopic 
Tefif-values in eq. 1 instead of the color temperatures. This 
would render the gravities larger by ~ 0.18 dex on aver- 
age, which is still lower than their spectroscopically de- 
rived values. It is feasible that this is due to a mild de- 
generacy of the entire parameter set, where a fraction of 
the offset can be explained by discrepancies in the adopted 
temperature. Since an increase in Tefi decreases the dif- 
ference between e(Fei) and e(Feii), one inevitably derives 
an increased stellar gravity. Due to this interdependence, 
a change in log 5 entails a slight change in the slope of 
the EP plot, which again is reflected by another rise in 
effective temperature (e.g., McWilliam et al. 1995; John- 
son 2002; Johnson et al. 2006: Fidbright et al. 2006; see 
also Section 3.3). Microturbulence also has an influence 
on log 5, though to a lesser extent. While matching e(Fei) 
and e(Fe 11) to achieve ionization equilibrium itself can be 
done accurately, one should also keep in mind that there is 
only a small number of stronger and well determined Fe II 
transitions available in our spectra and in the literature. 
On average, we have only about four Fell lines at hand, 
which will not only lead to larger random errors on the 
derived [Feii/H] values, but also frustrate accurate grav- 
ity estimates. As we will show in the next Section, even a 
implausibly large variation of \ogg has an only marginal 
effect on the abundances, in particular since we primarily 
rely on the neutral iron lines to assign the stellar [Fe/H]. 

The entire parameter determination for ^, Tefi and \ogg 
was successively iterated to approach simultaneous conver- 
gence, adopting the Fe i abundance from the last iteration 
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as the input metaUicity [M/H] (which strictly accounts for 
the global metal abundance, e.g., Salaris et al. 1993) for 
the subsequent atmosphere calculation. The final atmo- 
spheric parameters for each star derived in this way are 
tabulated in Tabic 3. With this final optimized parameter 
set, the abundance of the ionized state Fell deviates by 
no more than 0.04 dex from the Fei abundance, with an 
average deviation of -^0. 01-0. 02 dex. 

3.3. Abundance errors 

In order to determine accurate error estimates in the 
individual abundance ratios, we investigated the effects of 
the random uncertainties in the atmospheric parameters. 
First of all, abundance ratios were obtained for ten differ- 
ent model atmospheres for two stars (LG04c_000777 and 
LG04d_006628). These two stars were chosen because of 
their significantly different stellar parameters. In each of 
the models, the stellar parameters were separately varied 
by their typical random 1 a uncertainty, as conservatively 
estimated in the previous Section. In Table 4 the differ- 
ences of the abundance ratios for each of these models are 
given, i.e., for Toff ± lOOK, ^zbClkms"!, logg±0.2 dex 
and [M/II]±0.2 dex, w.r.t the values obtained from the 
atmospheres computed with the best spectroscopic param- 
eter set. 

At this stage of the analysis our model atmospheres for 
the Carina stars incorporate opacity distribution functions 
assuming scaled-solar element abundances (odfnew). It 
is important, however, to keep in mind that the a-elements 
are important electron donors in the atmospheres of red 
giants. These metals supply a significant number of free 
electrons, which sustain the formation of H~ the dom- 
inant source of the continuum opacity in these stars. It 
cannot be assumed a priori that the Carina stars exhibit 
any kind of enhancement in these elements, but measure- 
ments in seven dSphs of the Local Group, also comprising 
five stars in Carina (which do not coincide with our sam- 
ple), indicate that many stars in these galaxies are in fact 
enhanced in a number of a-elements (Shetrone et al. 2001, 
S03; Sadakane et al. 2004; Venn et al. 2004; Geisler et al. 
2005). Therefore we additionally derived abundances un- 
der the assumption of an a-enhancement of -1-0.4 dex in 
the atmospheres, by calculating models with the stellar de- 
fault parameters, now using the a-enhanced opacity distri- 
butions, AODFNEW (sec, e.g., Fulbright ct al. 2006, 2007). 
The change in each elements' abundance ratio derived with 
this atmosphere from the standard atmosphere's value is 
given in Table 4 in the column labeled "a" . 

The model parameters are not independent but rather 
intricately correlated in several ways (McWilliam ct al. 
1995). This leads to covariances of the parameters' uncer- 
tainties, which can be particularly significant for the case 
of Tefif and log 5 caused by an explicit temperature-gravity 
dependence. A change in Tes can also introduce a slope 
in the abundance vs. EW relation, which is used to fix 

To quantify the covariances of the stellar parameters, we 
follow the prescription of Johnson (2002), which is, in turn, 
based on the formalism of McWilliam ct al. (1995). For 
the representative stars LG04c.000777 and LG04d.006628 
we start by drawing 20 random temperature values from a 
Gaussian with their spectroscopic T^s as mean and with 
a dispersion of 100 K. For each of these, logg was changed 



accordingly, to reinstall ionization equilibrium. The re- 
sulting covariance 

N 

0-Tiogg = < T >){loggi- < log g >) 

amounts to 40, with a correlation coefficient of 0.79, i.e., 
the mutual influences of temperature and gravity on our 
abundances are in fact highly correlated. Likewise, we 
computed 20 values for the microturbulence with a dis- 
persion of O.lkms"^ and readjusted the gravity of the 
atmosphere to regain ionization equilibrium. These pa- 
rameters are less correlated, with cr^iogg = —0.005 and a 
correlation coefficient of 0.44. It turns out that covariance 
terms including the atmospheres' metaUicity [M/H] are 
negligible so that we did not include the respective mixed 
terms in the analysis. 

Finally, we account for random errors in the EW mea- 
surements by computing the la scatter in the abundances, 
weighted by the square root of the number of absorption 
features, N, used in the derivation of e (e.g., Carretta 
2006). This statistical error, ctew/V^, (sec also Table 5) 
accounts in principle for the line-to-line scatter and in- 
put atomic parameters. For those elements, for which 
only one absorption feature was measurable in the spec- 
tra, we adopted a representative error on the abundance 
through repeatedly varying the EW in the abundance de- 
termination within its typical accuracy (of the order of 
5-10%). In general, for those elements with iV > 5, which 
is the case for our Ca-, Til- and Fc-lines, the total abun- 
dance error is dominated by the uncertainties in the at- 
mospheric parameters (Prochaska et al. 2000). The ma- 
jor error source in terms of atmospheric parameters is in 
this case the tcmpciraturc sensitivity. For elements with 
few lines, there is a progressively larger contribution of 
EW errors due to enhanced noise, unresolved blends or 
uncertainties in the continuum placement. Interestingly, 
neglecting a-cnhancemcnt docs not significantly alter any 
of the derived abundance ratios and only has a significant 
effect on the ionized spcxdcs Ti 11 and Fe 11, where the en- 
hancement tends to overestimate the abundances by ~0.1 
dex. 

The total uncertainty ctot of each elemental abundance 
ratio [X/H] due to random errors in the atmospheric pa- 
rameters (listed in Table 5) is then obtained by summing 
the uncorrelated contributions from each individual pa- 
rameter (Teff, ^, logg, [M/H], a, AEW) in quadrature, 
and adding the respective covariant contributions from 
Teff, £, and logg (see eq. 2 in Johnson 2002). Since vari- 
ations in the atmospheres will affect different elements to 
a different extent, one has to additionally consider their 
covariances in computing the abundance ratios [X/Fe]. 
Therefore we obtain CTp-^Pgj = (Tx-l-cr|g—2crxFe with the co- 
variance (JxFo according to eq. 6 from Johnson (2002). As 
a result, the typical error on our derived [Fe I/H] amounts 
to ^0.17 dex and ranges from ^0.05 to ^0.20 dex for the 
other elements, while the uncertainty on [0/Fe] can be as 
high 0.35 dex. 

Finally, we note that none of our abundances has been 
corrected for NLTE effects. Since our results will be com- 
pared below with similar LTE analyses (e.g., S03; Venn 
et al. 2004), this source of uncertainty will not be system- 
atic on the whole. Despite our enforcement of ionization 
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equilibrium by reconciling e(Fe I) and e(Fe II), the neutral 
and ionized species of Ti show a marginal difference of 
—0.12 ± 0.18 dex. If we demanded ionization equilibrium 
for the Ti lines to establish a logg, our spectroscopic sur- 
face gravities would have had to be lower by 0. 13 dex on av- 
erage, which is well within our estimate of the uncertainty 
on this atmospheric parameter (see Section 3.2). The fact 
that our data do not exhibit any significant trend of £(Ti l) 
or £(Ti ll) with TeS nor with metallicity, supports our view 
that we can, for the moment, neglect NLTE effects in our 
analyses (Luck & Bond 1985; cf. Johnson 2002). 

4. RESULTS 

Table 5 shows the final abundance ratios [X/Fe] ([Fe/H] 
for the iron lines themselves) for each chemical element X. 
These ratios are given w.r.t. the mean abundance of the 
neutral iron lines, £(Fei). Also tabulated are the num- 
ber of lines, N, used to determine the unweighted mean 
abundance ratio for each element and standard deviation 
of each e measured from these lines. 

4.1. Iron 

The mean iron abundance of the red giants in Carina 
obtained from our high-resolution spectra is —1.69 dex 
with a standard deviation of 0.51 dex. This is in excel- 
lent agreement with the mean of —1.72 (on the scale of 
Carretta & Gratton 1997) derived from the CaT measure- 
ments. Also S03 find a mean of —1.64 with a scatter of 0.2 
dex from their high-resolution analysis of five red giants, 
while Smecker-Hane et al. (1999) reported a slightly lower 
value of —1.99 dex from a small number of stars with CaT 
measurements. Prom this point of view, our targets can 
he considered as Carina members representing the typical 
metallicity range in this dSph. 

There are some mildly metal-rich targets, which reach 
[Fe/H] of —1.27 dex. Furthermore, we detect two metal- 
poor stars with [Fe/II]< —2.5 dex, which contribute to 
the metal-poor tail of Carina's MDF. These stars lie 1.5- 
2 a below the peak of the MDF. Nevertheless, this small 
number of metal-poor objects is not sufficient to resolve 
the apparent G-dwarf problem (or in fact "K-giant prob- 
lem", Koch et al. 2007) in dSphs (Smecker-Hane et al. 
1999; Pont et al. 2004; Tolstoy et al. 2004; Venn et al. 
2004; Helmi et al. 2006; Paper I; Bosler et al. 2007; Koch 
et al. 2007). Interestingly, at 1.4 and 2 core radu, these 
two stars are also located at the largest radial distances 
in our sample, thus matching the mild radial abundance 
gradient that was detected in Paper I (sec also Harbcck et 
al. 2001). Recall that our parent sample selected for spec- 
troscopy was three times larger than the eventual member 
sample, and was selected to ensure that any very metal- 
poor or metal-rich member star would be included. The 
absence of such stars from our sample then indicates a real 
scarcity. 

4.1.1. Assessing the CaT calibration 

The six stars that were observed in both low- and high 
resolution mode enable us to assess the accuracy of the 
calibration of the CaT abundance scale onto iron, since 
these stars cover a range of abundances. Fig. 1 shows a 
comparison of the estimates from both methods. 



As the dashed line of unity in Fig. 1 (top panel) im- 
plies, all of the five more metal-rich red giants show good 
consistency between their CaT metallicity (on the scale 
of Carretta & Gratton 1997) and the high-resolution iron 
abundance, not deviating by more than 0.8 ct from unity. 
This range of metallicities is well sampled by the Galactic 
globular clusters (GGCs) in the calibrations of Rutledge et 
al. (1997a,b). An oft-cited caveat against calibrating the 
CaT in dSphs to a Galactic reference system is that GGCs 
tend to exhibit enhanced [Ca/Fe] ratios of ^ +0.4 dex on 
average (e.g., Pritzl et al. 2005), while dSph stars show 
lower values in this ratio (Venn et al. 2004). As we will 
show in Sect. 4.3, our Carina dSph stars exhibit slightly 
lower [Ca/Fe] ratios of the order of +0.09±0.14 dex, which 
should result in systematically lower line strengths of the 
CaT compared to the GGCs. A number of clusters in 
the Rutledge et al. (1997b) sample show in fact lower 
[Ca/Fe] values of ^ +0.1 dex at metallicities similar to 
that of our Carina stars (e.g., Cohen & Melendez 2005), 
and the GGCs, from which we derived the calibration of 
the CaT to [Fe/H]caT in Paper I, possess an average, en- 
hanced, [Ca/Fe] of +0.21+0.11 dex. Thus we might expect 
an inferred [Fe/H]caT lower by ^ 0.12 dex with respect to 
the [Fe/H] abundance, based on the difference in [Ca/Fe] 
between target and calibrators (see also discussions in Pa- 
per I; Bosler et al. 2007). We find that our CaT based 
measurements are more metal poor by (0.08+0.09) dex on 
average. This small zero point shift is consistent with the 
above difference in [Ca/Fe] of the calibration clusters and 
target stars. Note that there is no obvious trend of the 
difference between high-resolution [Fe/H] and [Fe/H]caT 
with [Ca/Fe] discernible in our data. 

The iron abundance of our metal-poor star is, at —2.50 
dex, lower by 0.44 dex (2.1 a) than the CaT value. This de- 
viation is in qualitative agreement with the CaT vs. high- 
resolution [Fe/H] comparison of a large number of dSph 
stars of Battaglia et al. (2008), who show that a system- 
atic overestimate of the CaT abundances sets in below 
[Fe/H]< —2.2 dex. Such a trend points to a general short- 
coming of the method to reliably infer iron metallicities 
from the CaT index in the most metal poor red giants. 
Since the globular clusters in the sample of Rutledge et al. 
(1997b) "only" comprised a metallicity range from —0.64 
down to —2.02 dex (Carretta & Gratton [1997] scale), mea- 
surements of stars beyond these limits necessarily rely on 
an extrapolation of the calibrations. In particular, the 
clusters employed for our calibrations in Paper I span a 
range from —2.0 to —1.12 dex. 

As an additional test, we calibrated our GGC measure- 
ments from Paper I onto the Fe II scale of Kraft & Ivans 
(2003). The major difference to the old calibration lies 
in the [Fe/H] value for the cluster M68, for which Kraft 
& Ivans (2003) find a more metal poor [Fe/H] of -2.43 
dex, while this cluster has a metallicity of —2 dex on 
the previously used scale of Carretta & Gratton (1997). 
In general, Carretta & Gratton (1997) have reported the 
most metal poor GGCs to be considerably more metal 
rich than other studies. The main impact of the new 
calibration (see bottom panel of Fig. 1) is that the most 
metal poor of our targets becomes even more deficient, at 
[Fe/H]caT ^ —3.1 dex, which underestimates the metallic- 
ity compared to our high-resolution measurement by 0.56 
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dex. On average, the CaT measurements on the Kraft & 
Ivans (2003) scale are more metal poor by 0.3 dex than the 
UVES data, with a standard deviation of 0.12 dex around 
unity. Although Cole et al. (2004) and, more recently, 
Carrera et al. (2007), have demonstrated that the CaT 
technique is generally applicable for a broad age range of 
0.25-13 Gyr and across a wide metallicity range spanning 
-2.2 <[Fe/H]< +0.47, the MDFs' mctal-poor tails remain 
unsatisfactorily sampled in most high resolution studies. 
The question of whether the metal-poor extensions of the 
MDFs of dwarf galaxies based on the CaT below ^ —2.2 
to —2.5 dex, which are reported throughout the literature 
(e.g., Smecker-Hane et al. 1999; Tolstoy et al. 2001; Bosler 
et al. 2007; Pont et al. 2004; Tolstoy et al. 2004; Cole et 
al. 2004, 2005; Paper I) are then reliable has to await more 
and highly accurate data in this regime (e.g., Battaglia et 
al. 2008). Nonetheless, the CaT index is well suited to 
the identification of candidates for fairly metal-poor stars, 
although the exact [Fe/H] value may not be measurable - 
this is then left for follow-up high-resolution studies. 

4.2. Alpha Elements - O, Mg, Si, Ca, Ti 

While iron is synthesized both in SNe II and SNe la, the 
even-Z elements (O, Mg, Si, Ca, Ti) are predicted to be 
produced during (nucleosynthetic) shell-burning in high- 
mass stars with very short lifetimes, which later end their 
lives in SNII explosions. 

A simple scenario for the formation of the Galactic halo 
is that an initial burst of SF produces large amounts of 
a-elements on short time scales, leading to an enhanced 
[a/Fe] ratio. At these early epochs, the long-lived SNe la 
did not yet have time to significantly pollute the initial 
SNe II abundance pattern. Once SNe la start to contribute 
(after a delay of ~1 Gyr from the onset of SF), iron pro- 
duction will increase, leading to a decline in the [a/Fe] 
abundance ratio. As noted above, the detailed patterns 
of the transitions in [a/Fe] with metallicity (taken as a 
proxy for time) are dependent on the specific SF history of 
a stellar system, on its stellar IMF, and on the timescales 
for mixing the SNe products into the interstellar medium 
(ISM) (Matteucci 2003; S03; Venn et al. 2004). 

The metal-poor stars of the Galactic halo evidently pop- 
ulate a plateau of [a/Fe]'^0.4 dex (see, e.g., the review by 
McWilliam 1997). From the available analyses of 36 red 
giants in seven dwarf galaxies, S03 found at best little ev- 
idence to identify these stellar populations with those in 
the Milky Way. In fact, stars in dSphs have been shown 
to be quite different from Galactic disk, and even more so 
from halo stars, in that they have lower [a/Fe] ratios (e.g., 
Shetrone et al. 2001; Pritzl et al. 2005; see also Fig. 3). 
Hence, these elements have often been considered as an 
important test for Searle & Zinn's (1978) hypothesis of the 
origin of the Galactic halo from satellite accretion events 
(Unavane, Wyse and Gilmore 1996). 

Despite the characteristic relatively low [a/Fe] in most 
stars of the dSphs, their wide abundance range from ap- 
proximately —3.0 to 0.3 dex and the fact that some stars 
also exhibit higher. Galactic [a/Fe] values suggests that 
using these elements as a test for chemical imprints of ac- 
cretion events requires a sophisticated analysis (Venn et al. 

The stellar abundances for these illustrations are tabulated in Venn 
sources. See their paper for references. 



2004). Irrespective of this analysis, the a elements provide 
important information in elucidating the evolutionary his- 
tories of individual dwarf galaxies (e.g., S03; T03). We will 
turn to discussing the implications of our derived [a/Fe] 
ratios for Carina's chemical evolution in Section 5. 

Figs. 2 through 4 display the derived abundances for 
each a-element in our Carina stars, together with the 
measurements of five red giants in Carina of S03, where 
none of their targets coincides with our present sam- 
ple. The a-elements are displayed individually as well 
as a combination of these elements as the parameter 
[a/Fe] = [(Mg-f-Ca-hTi)/3Fe] (see also Table 1 in Venn et 
al. 2004). The right panels of Fig. 2 and Figs. 3, 4 also 
include a compilation of Galactic thin and thick disk stars 
as well as a large number of halo field star^. Note that 
the separation into Galactic components is purely based 
on kinematic information. 

In order to compare Carina's elemental abundances with 
those of stars in other Local Group galaxies, we show 
also currently available abundance ratios for stars in eight 
Galactic dSphs (Draco, Sextans, UrsaMinor from Shetrone 
et al. 2001; Sculptor, Fornax, Leo I, Carina from S03; Ur- 
saMinor from Sadakane et al. 2004; Sculptor from Geisler 
et al. 2005 and Sagittarius from Monaco et al. 2005). 

When comparing data sets from various sources, as was 
done by Venn et al. (2004) and as is shown in Figs. 3 
and 4 here, one has to keep in mind that each analysis 
may use different spectral lines, presumably differing os- 
cillator strengths, and different model atmosphere tech- 
niques. The combination of such data will then result 
in an analysis-dependent spread in the abundance ratios, 
which is of the order of 0.1-0.2 dex (Venn et al. 2004), a 
value comparable to the abundance offsets found by S03 
by using different stellar atmospheres and atomic data. In 
the light of our measurement uncertainties (Sect. 3.3; indi- 
cated by the error bars in the figures) , we did not attempt 
to homogenize the data shown in the subsequent Figures. 
This measurement method induced scatter does not alter 
any of the conclusions we draw below. 

The first thing to note is that one of our target stars, 
LG04a_001826, exhibits a strong depletion in almost all 
of the a-elements, although it is moderately metal rich, 
at [Fe/H]=— 1.49 dex. This star will be discussed further 
below. 

As Venn et al. (2004) note, the scatter in [Mg/Fe] is ex- 
pected to be much larger than that in either [Ca/Fe] or 
[Ti/Fe], an expectation which is met by our data. Apart 
from the strongly discrepant target at [Fe/H]=— 1.49, with 
[Mg/Fe]=— 0.92, the more metal rich targets cover a full 
range in [Mg/Fe] from -0.24 to -1-0.08 dex, while the two 
most metal-poor stars in our sample exhibit enhanced Mg 
abundances of 0.3 and 0.4 dex, and thus partially over- 
lap with the Galactic halo. In order to test the extent to 
which the scatter can be considered intrinsic, we compared 
the observed scatter in the abundance ratios to the value 
predicted from a set of Monte Carlo simulations, in which 
our measurements were varied by the observational uncer- 
tainties. While the predicted and observed Ic-scatter are 
similar for the cases of [Ca/Fe] and [Ti/Fe], the observed 
scatter in [Mg/Fe], at 0.16 dex, exceeds that expected from 

al. (2004), who in turn assembled the data from various literature 
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the pure observational uncertainties by ~0.05 dex, indicat- 
ing a real dispersion. 

Ca is predominately produced in intermediate-mass pro- 
genitor SNell (Woosley & Weaver 1995). While the most 
mctal-poor stars of the Shctronc ct al. (2001) and 803 sam- 
ples in Sextans and Draco are clearly depleted w.r.t. the 
Galactic components (at [Ca/FeJ^O.l), our most metal- 
poor Carina giant is significantly enhanced in Ca (by 0.7 
dex). However, already the second most metal-poor star 
(0.2 dex more metal rich) has a [Ca/Fe] that is only slightly 
larger than the average we find fi-om the rest of our sample 
(< [Ca/Fe] > = 0.09±0. 14 dex). While the majority of our 
targets with [Fe/H]> —2 tend to follow the trend of the 
other dSphs' Ca abundances, and also are consistent with 
SOS's measurements, one should note that three or four of 
the more metal-rich giants have comparable enhancements 
to those of the Galactic thick disk and/or halo. 

Although it is still unclear if the nucleosynthesis of Ti 
proceeds in analogy to the other a-elements (e.g., Timmcs 
et al. 1995; Woosley & Weaver 1995), it is traditionally in- 
cluded in the group of these elements due to its similarly 
enhanced abundances in metal-poor stars (Gratton & Sne- 
den 1991). With the exception of the strongly depleted 
LG04a_001826, Ti in Carina tends to follow the overall 
trends of the other a-elements, with an average abundance 
significantly lower than the Galactic halo star abundance 
ratios. In fact, the two most metal-poor stars of our sam- 
ple show enhanced [Ti/Fe] of ^0.1 and 0.36 dex, which 
is consistent with the Galactic halo, but at odds with the 
two comparably metal-poor red giants from Shetrone et 
al. (2001) and S03. The latter stars, in the Sextans and 
Draco dSphs, are deficient in [Ti/Fe] by ~ —0.23 dex. The 
majority of our targets is consistent with both the Carina 
data and most of the other dSphs' abundances of Shetrone 
et al. (2001) and SOS. 

Si is often considered a prototypical Q-element. Its pro- 
duction site is identified with SNell of moderate progen- 
itor masses, around 20 M© (Woosley & Weaver 1995). In 
our data we were able to determine its abundance from 
typically three lines. For the most metal-poor star no in- 
formation could be extracted, since no measurable line was 
detected in its spectrum. 

In the second most metal-poor star, there is a strong in- 
dication of a significant enrichment in Si, as found in most 
metal-poor stellar abundance studies (McWilliam 1997; 
Prochaska et al. 2000; Venn et al. 2004, and references 
therein). As Fig. 4 implies, an enhancement of ~0.9 dex 
is not atypical for metal-poor giants in dSphs and also 
is found in UMi (Shetrone et al. 2001). Our stars show 
a lower average [Si/Fe] compared to the Shetrone sam- 
ples, where our mean value is <[Si/Fe]> = — 0.01 ± 0.12 
dex (for [Fe/H] > —2). It is worth noticing, however, that 
there is a fair amount of scatter in [Si/Fe], both in the 
Galactic halo and all dSph samples considered here, in- 
cluding our own data. In this context, [Si/Fe] abundance 
ratios as low as —0.4 are also found in other dSphs of the 
Local Group. While the second most metal poor star is 
enhanced in [Si/Fe], the abundance ratios around Carina's 
mean [Fe/H] typically reach from about —0.2 up to +0.15 
dex in our targets. Owing to the overall scatter, some of 
the dSphs' stellar abundances overlap with the Galactic 
halo component at [Si/Fe]~ +0.4. 



Oxygen abundances were measured from the two for- 
bidden lines at AA6300, 6363 A, which are, however, not 
detectable in the most metal poor target, nor in the overall 
depleted star of LG04a_001826. We find a range in [0/Fe] 
from —0.2 to +0.65 dex with a considerable star-to-star 
scatter at a given [Fe/H] (Fig. 5, top panels). As Fig. 5 
(bottom) implies, the range in [0/Fe] in Carina, as well as 
the increased scatter compared to the other a-elements, is 
consistent with that found in other dSphs. As mentioned 
above, these abundances are depleted w.r.t. the Galactic 
halo by ~0.5 dex and show a well defined decline with 
[Fe/H] at lower metallicities than in the Galaxy. 

4.3. Sodium 

Sodium abundances were determined from the absorp- 
tion lines at 5682 and 5688A. We did not include hyperfine 
splitting in the calculations, since its effect is expected to 
be neghgible (McWilliam et al. 1995). As Fig. 5 (bottom 
panels) implies, our mean [Na/Fe] in Carina of —0.36 dex 
with its scatter of ~0.27 dex is consistent with the Ca- 
rina values of SOS, and also with the [Na/Fe] in the other 
dSphs. In particular, there is a strong overlap of Carina's 
Na abundance with the Galactic halo stars, which also 
cover a broad rang(\ 

It is worth noticing that there is no trend of [Na/Fe] 
with oxygen discernible. Such a pattern would be evi- 
dence of deep mixing. In deep mixing, processed material 
from the interior of a star can be transported into the up- 
per atmosphere. Nucleosynthetic burning converts O, N 
and Ne to Na, whereas Mg is converted into Al in the 
deeper atmospheric layers of red giants. Hence, typical 
deep mixing patterns predict that an atmospheric deple- 
tion in O coincides with an enhancement of Na (Langer et 
al. 1993); while any observed Mg depletion is accompa- 
nied by an enhanced Al abundance (Langer & Hoffmann 
1995). Unfortunately, in none of our spectra could Al be 
detected, since all suitable lines in our spectral range are 
too weak. Interestingly, such (anti-) correlations are de- 
tected in globular clusters both in the Galaxy (Pritzl et 
al. 2005 and references therein) and in dSphs (Letarte et 
al. 2006), but they are not observed in Galactic field stars 
(Gratton et al. 2004) nor in dSph field stars (SOS; Venn 
et al. 2004). Consistent with this, our Carina red giants 
show a constant mean [Na/O] of ~ —0.63 ± 0.17 over the 
full range of metallicities. 

4.4. The a-depleted case of LG04a.001826 

We now turn to the question of whether our target star 
LG04a_001826, which shows a considerable depletion in al- 
most all of the a-clcmcnt ratios, i.e., in [Mg/Fe], [Ca/Fe] 
and [Ti/Fe], can provide further insight into general galac- 
tic evolution or whether it should rather be considered a 
special case. For this star [Mg/Fe] lies 2.2a below the sam- 
ple average, and [Ca/Fe] and [Ti/Fe] are underabundant 
by 1.8(7 and 2.3(t, respectively, w.r.t. the mean abundance 
ratios in Carina. 

One possible interpretation of a selective depletion of 
the a-elements could be the occurrence of deep mixing, 
but as discussed above, this is not seen in other stars in 
dSphs. Moreover, at [Na/Fe]=-0.48 LG04a_001826 shows 
a sodium abundance that is consistent with the mean 
[Na/Fe] ratio found from all our Carina stars, and Na 
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is not overly enhanced w.r.t. the non-detection of oxy- 
gen in its photosphere. It appears unhkely that this star's 
abundance patterns are influenced by deep mixing. An 
untestable hypothesis is that enrichment by rapidly ro- 
tating massive stars may be involved, since this can lead 
to globular-cluster-like abundance anomalies, dredging-up 
processed material into the outer layers of the massive 
stars, where it may be lost as a wind prior to the Type 
II supernova (Decressin et al. 2007). 

Generally speaking, the a-clcmcnts arc synthesized in 
core-collapse SNell, where difl'erent elements originate in 
different stages of these endphases in the lives of massive 
stars. Oxygen and magnesium, which show the strongest 
depletion (Mg) or upper limits (O) in this star, are 
produced mainly during the hydrostatic burning phases. 
Their yields are primarily a function of the progenitor's 
mass, but are heavily affected by the thickness of the burn- 
ing shells, the treatment of rotation and convection, and 
the mass cut between the proto-neutron star and the SN 
ejecta (Thielemann et al. 1996; Argast et al. 2002). If one 
assumes a larger range of possible realizations of these fac- 
tors, a larger scatter in the resulting element ratios such as 
[0/Fe] and [Mg/Fc] is feasible (Argast et al. 2002; Venn et 
al. 2004). Depending on the exact mass of the SNII pro- 
genitor, a low value of [Mg/Fe] is possible (see Fig. 1 in Ar- 
gast et al. 2002). Low [Mg/Fe] is observed in a few metal 
poor bulge-stars (McWilliam et al. 1995) and in mildly 
metal-poor halo stars (e.g., Carney et al. 1997; Ivans et al. 
2003). 

Some SN models (like those of Argast et al. 2000, 2002) 

predict the occurrence of such depleted stars at lower, 
halo-like metallicities, specifically due to the mixing of 
SN ejecta into a gas reservoir of ^ lO^M©. Compara- 
ble variations at a higher [Fe/H], as seen in the case of 
LG04a_001826, can be explained analogously with reason- 
able parameters, particularly noting the uncertainties as- 
sociated with the modelling of Mg and Fe yields. Alter- 
natively, one could assume that our target star exhibits a 
normal average abundance of Mg, but has an anomalous 
Fe abundance. Then its low [Mg/H] would be represen- 
tative of a halo-like star of [Fc/H]^ —2.6 dex. In this 
case, in order to realize its comparatively low [Mg/Fe], 
LG04a_001826 had to originate from SNela pre-enriched 
material (Gilmore & Wyse 1998). This would be diluted 
less than the characteristic halo ISM or could be subject to 
a significantly higher iron yield, producing a [Fe/H] higher 
by ~ 1.1 dex than expected from its a-yield. 

The other a-elements Si, Ca and Ti originate in explo- 
sive burning during the SN event itself and depend primar- 
ily on the explosion energy and to a lesser extent on the 
amount of material available to be burned and the bind- 
ing energy of the star, that is, its mass (e.g., Woosley & 
Weaver 1995). Since typical explosion energies cover only 
a small range, the abundance ratios of these elements can 
be expected to show a much lesser scatter, which is in fact 
in agreement with the observations (Venn et al. 2004; see 
also Figs. 3 and 4). 

It is then odd that this star exhibits a rather average 
[Si/Fe] ratio (at 0.2(t below the sample mean). It has been 
suggested that Si contains a non-negligible admixture of 
material originating in hydrostatic burning (Timmes et al. 
1995; Thielemann et al. 1996; Argast et al. 2000) and it 



is then likely that the resulting abundance patterns can 
differ from a pure explosive origin as for Ca or Ti (but see 
Cohen et al. 2007 for the exact opposite pattern). 

The intermediate-mass a-element Ti has also a contri- 
bution from complete Si burning with an a-rich freeze-out 
(Thielemann et al. 1996 and references therein). Differ- 
ences between Ca, Si, and Ti are thus naturally to be 
expected. Finally, the production of Ti is not yet fully 
understood (Woosley 1986; Woosley & Weaver 1995) and 
Galactic evolution models often fail to reproduce the ob- 
served abundance patterns (e.g., Timmes et al. 1995). A 
detailed interpretation of this star's abundance ratios in 
terms of the yields for each individual a-clcmcnt has to rely 
on fine-tuned nucleosynthetic modeling and is presently 
hampered by the lack of an accurate understanding of the 
exact SNII explosion mechanism (e.g., see discussions in 
Woosley & Weaver 1995; Timmes et al. 1995; McWilUam 
et al. 1995; Thielemann et al. 1996; Argast et al. 2000, 
2002). 

All in all, the presence of a strongly a-element depleted 
object would seem to suggest that we see here the influ- 
ence of single SN events. Contrary to a homogeneous, 
close-to-uniform [a/Fe] pattern resulting from well- mixed, 
IMF-averaged SNe, this star is likely to reflect the imprint 
of stochastic inhomogeneous chemical evolution, where the 
star formed from unmixed SN ejecta. Rather than (over- 
) interpreting this star's abundance pattern in terms of 
Carina's SF history we attribute it to statistical fluctu- 
ations due to incomplete ISM mixing. That such incom- 
plete mixing may occur also in actively star-forming dwarf 
galaxies was recently shown by Kniazev et al. (2005). It 
is then also possible to assign the strong Ca enhancement 
of the most metal-poor star in our sample to such scat- 
ter, which becomes progressively dominant in metal-poor 
stars of [Fe/H]< —3 dex. This contrasts with the low scat- 
ter in the more metal poor Galactic halo, indicative of a 
well-mixed ISM even at very early times (e.g., Cayrel et 
al. 2004; Frangois et al. 2004). Of course, any distinction 
between Carina and the halo based on the small number 
of peculiar stars in Carina should be taken with caution: 
while the dispersion in the abundance ratios of metal poor 
halo stars is in fact smaller than would be expected if all 
metal poor stars were polluted by single SNe, there are a 
number of chemically peculiar metal poor stars present in 
the halo (e.g., Norris et al. 2002; Cohen et al. 2007). 

Fortunately, observational analyses of a-deplctcd stars 
in the Galactic halo are underway, and will quantify how 
SNela and SNe II contribute to their abundance pat- 
terns. For instance, Ivans et al. (2003; and references 
therein) investigate three moderately metal poor halo stars 
([Fe/H]^ —2) with unusually low a- and neutron capture 
element abundances. By comparison with the abundance 
patterns of very metal poor stars, these stars were argued 
to be iron rich rather than a-poor. As in our Carina case, 
their depletion in [a/Fe] would be explicable in terms of a 
contribution of SNe la yields that is larger than expected 
given their low metallicities. It is conceivable that such 
stars may have been born out of material that was pol- 
luted by the earliest SNe la to have occurred in the Milky 
Way (Nissen & Schuster 1997; Ivans et al. 2003). It has 
been suggested that at least a fraction of this new class 
of Galactic halo stars may have been captured from ac- 
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creted, dSph-like systems. Despite the general discrepancy 
between the halo and present-day dSph abundances, our 
finding of such a peculiar, a-depleted star in Carina then 
makes feasible a partial accretion origin of this chemically 
distinct halo subpopulation. 

5. IMPLICATIONS FOR CARINA'S CHEMICAL EVOLUTION 

Apart from the two most metal-poor stars in our sam- 
ple and the strongly depleted red giant LG04a_001826, 
the majority of our targets are compatible with solar a- 
abundances: we find a mean [a/Fe] of O.OOiO.lO dex. 
Our data confirm the earlier important finding that the 
a-elements in dSph stars are significantly depleted w.r.t. 
Galactic halo stars of the same mctallicity (Shctronc et al. 
2001; S03; T03; Venn et al. 2004). It is interesting to note 
that the two most metal-poor stars in Carina show an en- 
hancement in the a-elements by ^ 0.3 and 0.4 dex, which 
is the same order of magnitude as seen in the halo stars. 
These stars should then also be expected to be the oldest 
stars to have formed and may have been born at an early 
epoch, when SNe la had not yet started to contribute, but 
when an initial starburst had already polluted the pris- 
tine material with a significant amount of the a-elements 
ejected in short-lived SNell. In general, dSph stars pop- 
ulating the same SNe II plateau as the Galactic halo are 
consistent with an invariant, massive star, IMF (e.g., Wyse 
ct al. 2002). 

On the other hand, those stars with higher mctallicities 
(between ~ —1.8 and —1.3 dex) resemble stars in other 
dSph galaxies. These progressively lower values of [a/Fe] 
are characteristic of the slow chemical evolution and low 
SF rates in dwarf galaxies (Unavanc ct al. 1996; Matteucci 
2003). Differences between individual galaxies' abundance 
ratios are, of course, to be expected, since these chemical 
tracers sensitively depend on the details of each galaxy's 
environment, SF history, IMF, etc. Overall, our stars' 
abundance ratios are also in good agreement with the val- 
ues found in the five stars studied by SOS. 

It has been proposed that Carina's episodic SF history, 
i.e., its sequence of bursty and quiescent phases, is re- 
flected in the [a/Fe] abundance ratios (S03, T03). 803 
suggested that the considerable rise in [a/Fe] - after the 
typical decline due to the continuing feedback from SNe la 
- seen in their data at [Fe/H]~ —1.65 indicates that these 
stars may be related to cither Carina's second SF episode 
around 7 Gyr ago or to an intermediate SF event well be- 
fore this episode. Both interpretations are compatible with 
photometric determinations of Carina's SF history and the 
times of each SF episode and hiatus (Smecker-Hane et al. 
1996; Hurley-Keller et al. 1998; Rizzi et al. 2003). Any 
connection is strengthened by comparable variations in a 
number of heavy clement abundance ratios. Resolving and 
quantifying the SF history clearly requires further inves- 
tigation: the interpretation of S03 and T03 was based on 
only 5 stars with limited age resolution. 

We combine our data for the a-elements O, Mg, Si and 
Ca with the measurements from S03. The resulting data 
set represents the entire currently available high-resolution 
abundance data for Carina (we do have a larger study un- 
derway). It is tempting to interpret the [a/Fe] pattern 
of the full data set (Figs. 3 and 6) in terms of Carina's 
episodic SF. S03 have noted a rise in [a/Fe] between —1.65 



and —1.6 dex from their 5 stars. This rise is also dis- 
cernible in the full data set including S03's and our ten 
stars, where it is most distinct in [0/Fe] and less pro- 
nounced in the other elements. After a decline in SF rates, 
the [a/Fe] ratio may rise again towards a maximum at 
around [Fe/H]=— 1.4 dex, before levehng off again. How- 
ever, in the light of the measurement uncertainties (see 
also Fig. 6, bottom right) the run of [a/Fe] vs. [Fe/H] for 
most of our stars (disregarding the most metal-poor and 
strongly a-dcplctcd star for the moment) and those of S03 
can equally well be considered to be constant (see bottom 
right panel of Fig. 6), with some, perhaps intrinsic, scatter 
in [a/Fe]. Uniformly low [a/Fe] ratios over a wide range 
of mctallicities would then point to the fact that Carina's 
SF was in fact dominated by SNe II driven blow-out (e.g., 
Mac Low & Ferrara 1999; Paper I). 

In Fig. 7, we compare the full observational set with the 
models from Lanfranchi ct al. (2006, hereafter LMC06), 
which are an extension of the models presented in Lan- 
franchi & Matteucci (2004). This set of models was tai- 
lored, modifying best-fit SF and wind efficiencies, so as 
to reproduce the total mass, present-day gas mass and the 
few elemental abundances available at that time for a num- 
ber of dSphs. In addition to the aforementioned observ- 
ables used in Lanfranchi & Matteucci (2004), LMC06 also 
incorporated our MDF of Carina from Paper I to update 
their choice of model parameters. A marginally higher SF 
efficiency (of 0.15 Gyr^^) and a slightly lower wind effi- 
ciency (of 5 times the SF rate) were adopted to better 
match our observations. The models described in Lan- 
franchi & Matteucci (2004) assumed Carina's SF history 
according to Hernandez, Gilmore, & Valls-Gabaud (2000), 
and so are lacking a significant very old population, (e.g., 
Smecker-Hane et al. 1994; Harbeck et al. 2001; Monelli et 
al. 2003; Grebel & Gallagher 2004). The LMC06 models 
adopt the SF history from Rizzi et al. (2003). In this lat- 
ter study, the SF was found to have occurred episodically 
with an initial active period at an early epoch, while the 
majority of stars were formed ~ 3-7 Gyr ago. 

For all the elements shown in Fig. 7 there is a sudden 
decrease in the computed [a/Fe] at [Fe/H]^ —1.6, consis- 
tent with a combination of the onset of galactic winds and 
the injection of Fe into the ISM via SNIa events. Once 
the winds set in, the available amount of gas for SF is 
reduced, resulting in a decline of the overall SF. Conse- 
quently, the formation of new stars almost ceases and, in 
turn, the occurrence of SNe II and the related production 
and injection of the a-elements is halted as well. Since 
Fe is continuously produced and injected into the ISM by 
the longer timescale SNe la for several Gyr, unaffected by 
cessation of SF, the start of efficient winds leads to the sud- 
den decrease seen in the modeled [a/Fe] abundance ratios. 
It is, however, not obvious a priori that the SF rate and 
thus the associated SNe II rate was high enough to sustain 
winds - these free parameters are put in by hand in the 
LMC06 models. 

For Ca the majority of our stars lie above the models, 
while the observed [Mg/Fe] and [Si/Fe] ratios tend to lie 
below the predictions at any given mctallicity. This may 
in part be attributed to the large scatter in [a/Fe] for 
basically all metallicities, a scatter which also has been re- 
ported to exist in previous observational studies (see Pritzl 
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et al. 2005 and references therein). At present, the models 
appear to provide at best an upper limit to the observa- 
tions, given the overall large star to star scatter, and thus 
reproduce the more chemically enhanced stars reasonably 
well. One should keep in mind that the original model pa- 
rameters were optimized to fit the abundance patterns as 
observed by S03 so that the overall trends in their element 
ratios are fit fairly well. 

The overall impression to glean from Fig. 7 is that, in 
spite of being able to predict certain patterns in chemi- 
cal abundance ratios, there is a general need for models 
of chemical evolution to steadily improve the parameters 
governing the SF processes. In essence, matching the full 
sample of available a-elemental data appears to require 
a turnoff from the plateau in [a/Fe] at lower metallici- 
ties (again with the exception of Ca). These low abun- 
dances at higher metallicities may be achieved by adopt- 
ing a less intense wind and/or a cessation in SF at earlier 
epochs. More intense winds would provide lower values 
for all [a/Fe] ratios, once the wind developed, due to the 
correspondingly decreased SF rate. Limitations can be 
attributed to the difficulties in assuming a specific SF his- 
tory for the models and will depend both on the number, 
duration, and times of the assigned SF episodes. Simi- 
larly, an overall lower SF efiiciency would be able to re- 
tain the lower [a/Fe] values. Since these parameters affect 
the shape of Carina's MDF (recall that they have been 
adjusted to match the latter) it is clear that ideally the 
MDF, the full set of observed abundance ratios and a fully 
sampled SF history have to be simultaneously taken into 
account in order to reliably devise the evolutionary his- 
tory of Carina. Considering the discussion in Sect. 4.4, it 
is also conceivable that stochastic SF, e.g., in small associ- 
ations, and inhomogeneous enrichment, may have played 
a major role in Carina's history. 

In Fig. 8 we overplot a set of chemical evolution models 
from Gilmore & Wyse (1991) with our observed oxygen 
abundance ratios. In these models, which were the first 
chemical evolution models to account naturally for below- 
solar a-element ratios at low metallicities, SF occurs in 
an initial, short-lived burst, during which there is no sig- 
nificant contribution from SNe la. This active period is 
followed by a quiescent phase lasting several Gyr, during 
which only SNe la enrich the ISM. A second, late SF burst 
is assumed to set in at higher metallicities. In the figure, 
the first burst occurs at [Fe/H] =— 2.1 dex, while the sec- 
ond SF sets in at —1.1 (solid line) and —1.4 dex (dashed 
line). We did not aim at optimizing the exact onset of each 
burst to constrain Carina's SF history, but rather shifted 
the predictions of Gilmore & Wyse (1991), which were 
originally computed to represent the chemical evolution of 
the LMC and SMC. It is then intriguing that already this 
simple scenario yields an amazingly good representation 
of our data, and that the observed increase and scatter 
in the abundance patterns can in fact be assigned to an 
increase in SF activity, perhaps associated with separate 
SF episodes 7 and 3 Gyr ago (Hurley-Keller et al. 1998). 
This comparison indicates that the onset of the second 
epoch of SF has to be around an abundance of —1.5 dex. 
It cannot be ruled out by our data that two SF bursts 
have occurred in the interval between ~ —2 and ~ —1.4 
dex, which is probable given the fact that deep photo- 



metric studies agree on the occurrence of three major SF 
episodes in Carina (Smecker-Hane et al. 1996). 

The occurrence of SF bursts should cause variations in 
both [0/Fe] and the other a-elements, whereas the ob- 
served scatter in [Ca/Fc] and [Ti/Fe] is smaller than that 
seen in our [0/Fe] abundance ratios. As Gilmore & Wyse 
(1991) argue, inhomogeneities in the location and onset 
of each SF burst inevitably leads to a large scatter in the 
abundance ratios, such as seen in Carina. Even in the pres- 
ence of a statistically universal IMF, the intrinsically rare 
high-mass stars will form with significant random number 
fluctuations in SF events of small total mass: SF events 
that produce only a small total mass in stars might not 
form any of the highest-mass oxygen-producing stars at 
aU. Adopting a standard IMF (e.g., Miller & Scalo 1979), 
it can be shown that Mtot ~ lOOM© form in low-mass 
stars per SNII (M > lOM©), while ~ 4000^4; form in the 
lower-mass regime of the IMF for every 4OM0 star. Thus, 
in a system like Carina, where SF involves small num- 
bers, as well as low rates, the resulting mass-averaged IMF 
yield will have inevitable stochastic scatter, contributing 
to observed intrinsic scatter in our abundance ratios. Ef- 
fective truncation of the massive-star IMF due to a low 
SF rate was suggested by T03 as a mechanism to pro- 
duce low values of [a/Fe] in dSph stars, treating the entire 
galaxy as a "single-cell star-forming entity". Our model 
comparison and the new, more numerous data allow us 
here to identify additional new key elements inhomoge- 
neous, stochastic SF in low-mass associations or clusters, 
with incomplete mixing leading to inhomogeneous chem- 
ical evolution. Overall, our data suggest that Carina's 
evolution involved non-monotonic, stochastic, low-rate SF 
that forms marginally bound low-mass stellar associations 
and/or clusters, convolved with incomplete mixing of the 
SNe ejecta into the ISM. 

6. CONCLUSIONS 

We have derived elemental abundance ratios for a num- 
ber of light chemical elements in ten red giants in the Ca- 
rina dSph, tripling the existing sample for this galaxy. By 
adding these ten stars to the presently existing data for 
~50 red giants in eight out of the 17 currently known dSph 
satellites of the MW system, the sparse chemical informa- 
tion available for dSph galaxies slowly continues to grow, 
confirming the overall trends in element abundance ratios 
noticed in earlier studies. 

As a number of our stars were also observed at lower 
resolution, we could compare metallicity estimates from 
the well established calcium triplet (CaT) calibration with 
the high-quality iron abundances from the present work. 
We show that the CaT-based metallicities and the true 
iron abundances agree remarkably well in the relatively 
metal-rich regime of [Fe/II]> —2 dex. For the single more 
metal-poor star overlapping with the high-resolution sam- 
ple wc find a substantial difference in the metallicity de- 
rived from the CaT and high-resolution iron lines. This 
may be since the CaT calibration needs to be extrapolated 
below [Fe/H]< —2 dex. Since we have only one data point, 
a more complete sampling of the dSphs' metal-poor stars 
with high-resolution spectroscopy is desirable. The red gi- 
ants in Carina span a comparable range in [Ca/Fe] to the 
GGCs that we used to calibrate Carina's low-resolution 
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MDF in Paper I, so we argue that the apphcation of the 
CaT method to devise a global metallicity is valid within 
the range of the calibration. 

The majority of the elements observed follow the pat- 
terns seen in other MW satellite dSphs (SOS; Venn et al. 
2004), i.e., they lie generally below the Galactic halo star 
values. This is consistent with low SF efficiency governing 
the internal evolution of the dwarfs, but is no consistent 
with halo formation by accretion of present-day satellites 
(Shetrone et al. 2001). Our Mg abundances show larger 
scatter than do Ca and Si, as anticipated by their different 
nucleosynthetic formation processes. 

We detect two metal-poor objects at —2.5 and —2.7 dex. 
These are significantly a-enhanced and exhibit halo-like 
abundance ratios, while the most metal-poor star is even 
more strongly enhanced in Ca. Overall, it appears that the 
Ca ratios are slightly higher on average than the remain- 
der of the a-elements. Therefore the relatively metal poor 
stars in the Galactic halo can in fact be considered as rem- 
nants from an accreted dSph satellite, if this merger event 
occurred very early on (e.g., Grebel et al. 2003; CaroUo et 
al. 2003). 

Our sample hosts one star that is remarkably depleted 
in almost all elements. Possible explanations include a 
particularly high iron yield or the imprint of single SN 
events, the latter of which points to the presence of early 
stochastical fluctuations and an incomplete mixing of the 
ISM gas. This may indicate that SF in Carina occurred 
stochastically within smaller constituents, followed by an 
inhomogeneous enrichment of the surrounding medium. If 
SF in dissolving star clusters was a major contributor in 
Carina, one might be able to identify the clusters' rem- 
nants. Chemical substructure is, however, generally not 
observed in the low mass dwarf galaxies (but see Kleyna 
et al. 2003), nor in transition-type galaxies. Localized ev- 
idence of abundance imhomogeneities on small scales at 
a given time, which argues in favor of SF in associations, 
has been detected in a few dwarf irregular galaxies (e.g.. 



Grebel & Richtler 1992; Gonzalez & Wallerstein 1999; Kni- 
azev et al. 2005). 

We compare our data to models that include several 
bursts in Carina's SF history and conclude that the an- 
ticipated low SF efficiencies and moderate mass-loss rates 
in dSphs provide reasonable representations of present-day 
data sets, although they fail to produce the observed abun- 
dance scatter. There are still deficiencies in the detailed 
models concerning some single elements. In particular the 
Ca measurements are at odds with the models, having 
systematically higher abundances on average. We note 
that the LMC06 models that we compared to our data 
did not account for incomplete mixing, which we, how- 
ever, observe. By adopting a model prediction with two 
SF bursts and an extended hiatus, which is qualitatively 
justified by Carina's SF history from its CMD, we argue 
that the second burst in Carina must have occurred at an 
approximate [Fe/H] of —1.5. This burst metallicity is not 
quantitatively constrained by CMDs, and hence is a pre- 
diction for future abundance measurements (see also S03). 
Such predictions become testable, once metallicities and 
robust age estimates of stars on the subgiant branches of 
Carina are available. Although our data do not rule out 
any imprints of Carina's episodic SF, we did not attempt 
to quantify the underlying exact SF history, which is fur- 
ther impeded by the lack of accurate age information for 
our target stars. 
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Fig. 1. — Comparison of the mctallicitics froiii high- and medium-resolution (CaT) analyses on the metallicity scales of Carretta & Gratton 
(1997; CG97; top panel) and Kraft k, Ivans (200.3: K103: bottom ijancl). Shown arc those six of our target stars, for which CaT widths were 
measured in Paper I. The horizontal crrorbars are formal measurements errors. The 1 u-errorbars on [Fe l/H] are those deduced in Section 3.3. 
The grey-shaded areas illustrate the ranges, where the CaT calibrations have to be extrapolated due to the lack of very metal poor calibrating 
clusters in Rutledge et al. (1997a,b) and Kraft & Ivans (2003). 
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Fig. 2. — Abundance ratios for the a-elements Mg, Ca and Ti: Our measurements for the ten Carina red giants are shown as filled blue 
stars, the filled black squares in the left panels are the Carina stars from 303. The additional data points in the right panels are Galactic field 
stars, where the color coding illustrates the kinematic separation into Galactic components and is shown in analogy to Venn et al. (2004); 
halo (cyan), thick disk (green), thin disk (red), retrograde orbits (black) and those with a high- velocity Toomre component (blue). For stars 
shown as open cyan symbols, there is no existing velocity information. The errorbars in the upper right corner are those derived in Section 3.3. 
[See the electronic edition of the Journal for a color version of this figure.] 
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Fig. 3. — Average a abundance ratio in analogy to Fig. 2. In addition, chemical element data of stars in eight other dSphs are shown as 
open squares (Shetrone et al. 2001, S03; Sadakane et al. 2004; Geisler et al. 2005; Monaco et al. 2005). [See the electronic edition of the 
Journal for a color version of this figure.] 
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Fig. 4. — Same as Fig. 3, but for silicon abundances. Due to the lack of Si information in Venn et al. (2004), this figure was taken from 
Pritzl et al. (2005) and additionally displays a sample of globular clusters (filled large circles). Reproduced with kind permission by B. J. 
Pritzl. [See the electronic edition of the Journal for a color version of this figure.] 
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Fig. 5. — Same as Fig. 2, but for Na and O abundances. For illustration, the oxygen data were complemented with halo stars from the 
samples of Israelian et al. (1998); Boesgaard et al. (1999); Mishenina et al. (2000); Akerman et al. (2004); Bai et al. (2004); Cayrel et al. 
(2004) and Zhang & Zhao (2005). [See the electronic edition of the Journal for a color version of this figure.] 
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Fig. 6. — Our measured [o/Fe] (blue stars) shown together with the data from SOS (black squares). The bottom panels are blow-ups of the 
region enclosed by the dashed quadrangle in the top panel. In the bottom left panel errorbars were omitted for clarity. [See the electronic 
edition of the Journal for a color version of this figure.] 




Fig. 7. — Our measured abundance ratios for the a-elements O, Mg, Si and Ca, shown as blue stars, in comparison with the model 
calculations from LMC06 (solid lines). Also plotted are the five data points from SOS (black squares). [See the electronic edition of the 
Journal for a color version of this figure.] 
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Fig. 8. — Chemical evolution models from Gilmore & Wyse (1991), which incorporate an initial SF burst, an extended quiescent phase, and 
a second successive burst. The solid and dashed line adopt different onsets of SF. The models were offset by —0.4 dex in [Fe/H] w.r.t. the 
original predictions from Gilmore & Wyse (1991) to match the observed element ratios from this work (blue stars) and SOS (black squares). 
[See the electronic edition of the Journal for a color version of this figure.] 
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Table 1 

Radial velocity members in Carina. 



Star" 


tt (.12000) 


() (.12000) 


\" 


B-V 


\-K 


\"r [klllb 


S/N 


LG04b.004260 


06 


40 


30 


-50 


59 


15 


17.75 


0.99 


3.17 


215.6 


18 


LG04C.006477 


06 


41 


21 


-51 


03 


43 


17.71 


1.31 


2.71 


219.5 


17 


LG04a.000377 


06 


41 


59 


-50 


51 


13 


17.78 


1.30 


3.23 


219.5 


18 


LG04a_001556 


06 


42 


17 


-50 


55 


55 


17.82 


1.25 


3.10 


223.1 


19 


LG04a_001826 


06 


40 


58 


-50 


53 


35 


17.76 


1.29 


2.88 


224.8 


18 


LG04C.000777 


06 


41 


04 


-51 


01 


35 


17.55 


1.31 


3.02 


225.1 


20 


LG04a_002181 


06 


41 


39 


-50 


49 


58 


17.62 


1.36 


2.88 


230.6 


21 


LG04d_006628 


06 


41 


37 


-51 


01 


43 


17.95 


1.28 


2.90 


231.1 


23 


LG04c_000951 


06 


39 


55 


-50 


57 


36 


17.79 


1.19 


2.72 


233.2 


15 


LG04C.000626 


06 


40 


47 


-51 


06 


03 


17.60 


1.37 


3.10 


237.6 


16 



Note. — The given colors are dereddened values in the JC system. 

"The nomenclature is such that LG04a-d designates the ElS-ficlds targeting Carina, followed 

by the number in the respective input catalog (P.Lynam [ElS-team], 2003, private communication). 



Table 2 

LiNELIST 



Ion 


Wavelength 

[A] 


EP 
[eV] 


log gf 










Equivalent Width 


[mA] 










Arcturus 


4260" 


6477" 


377" 


1556" 


1826" 


777" 


2181" 


6628" 


951" 


626" 


Pel 


4802.88 


3.69 


-1.64 


89.8 


38.8 


35.0 


79.3 


20.5 


37.5 


39.6 


43.5 


21.5 






Fel 


4839.54 


3.27 


-1.90 


102.4 


.34.0 


42.1 


69.8 


44.6 


60.2 


54.3 


59.5 


30.5 


29.3 


38.1 


Fel 


4871.32 


2.85 


-0.43 


246.9 


91.8 


102.5 


120.8 


126.6 


114.5 


125.6 


141.1 


82.7 


74.5 


116.2 


Pel 


4903.31 


2.88 


-0.67 


206.6 


88.9 


107.7 


109.9 


107.9 


95.6 


97.0 


112.4 


71.7 


36.1 


70.2 


Pel 


4918.99 


2.85 


-0.41 


250.5 


136.1 


127.2 


142.9 


137.5 


141.2 


151.0 


142.4 


95.0 


74.2 


98.4 



Note. — This Table is published in its entirety in the electronic edition of the Astronomical Journal. 

"For the saJje of readability, only the catalogue-numbers are given here, omitting the prefixes "LG04[abcd]_" (cf. Table 1) 



Table 3 

Spectroscopic and Photometric model atmosphere parameters 



Star 


Teff 




Teff 




log g 


IFe/H] 


[Af/H] 


Photometric 


Spectroscopic 


[km / s] 


cgs 


(CaT) 


LG04b_004260 


4091 


62 


4594 


1.52 


1.53 


-1.63 


-1.50 


LG04C.006477 


4397 


80 


4600 


1..54 


1.49 


-1.55 


-1.42 


LG04a_000377 


4060 


60 


4511 


1.74 


1..39 


-1.45 


-1..50 


LG04a_001556 


4134 


64 


4570 


1.75 


1.68 


-1.41 


-1.30 


LG04a_001826 


4271 


73 


4384 


1.74 


0.84 




-1.46 


LG04C.000777 


4183 


67 


4205 


1.25 


0.85 




-1.38 


LG04a_002181 


4279 


73 


4317 


1.66 


1.26 


-1.39 


-1.39 


LG04d_006628 


4260 


72 


4717 


1.40 


2.02 




-1.81 


LG04c_000951 


4387 


80 


4395 


1.66 


0.90 




-2.74 


LG04c_000626 


4130 


64 


4420 


1.70 


1.80 


-2.06 


-2.57 
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Table 4 

Error analysis for two representative stars - uncorrelated contributions 



STW mSid SI A [JW/HJ 

+100K -lOOK +0.2 -0.2 +0.1 -0.1 +0.2 -0.2 

LG04C.000777 



Oi 


-0.010 





+0.090 


-0.080 


-0.010 





+0.060 


-0.070 


+0.080 


Nai 


+0.090 


-0.095 


-0.035 


+0.040 


-0.010 


+0.015 


-0.025 


+0.045 


-0.045 


Mgl 


-0.039 


+0.028 


-0.018 


+0.017 


+0.049 


-0.051 


-0.012 


+0.011 


-0.036 


Si I 


-0.132 


+0.138 


+0.032 


-0.033 


+0.069 


-0.065 


+0.025 


-0.023 


+0.010 


Cai 


+0.023 


-0.054 


-0.025 


+0.023 


+0.021 


-0.021 


-0.035 


+0.027 


-0.037 


Til 


+0.094 


-0.125 


-0.006 


+0.011 


+0.033 


-0.030 


-0.034 


+0.029 


-0.035 


Tin 


-0.130 


+0.113 


+0.056 


-0.059 


+0.019 


-0.015 


+0.041 


-0.045 


+0.056 


Fel 


+0.116 


-0.081 


-0.002 


+0.003 


-0.075 


+0.075 


+0.005 


-0.001 


+0.040 


Fen 


-0.112 


+0.162 


+0.065 


-0.070 


-0.053 


+0.058 


+0.070 


-0.062 


+0.125 


LG04d_006628 


Oi 

Nal 


-0.010 


+0.010 


+0.080 


-0.090 








+0.070 


-0.070 


+0.090 


Mgl 


-0.040 


+0.054 


-0.011 


+0.023 


+0.010 


-0.003 


+0.006 


-0.006 


+0.012 


Si I 


-0.110 


+0.114 


+0.029 


-0.027 


+0.040 


-0.043 


+0.016 


-0.026 


+0.032 


Cal 


-0.044 


+0.039 


-0.001 


+0.003 


+0.020 


-0.020 


+0.006 


-0.009 


+0.001 


Til 


+0.022 


-0.041 


-0.005 


-0.003 


+0.024 


-0.025 


-0.001 


+0.005 


-0.015 


Ti 11 


-0.115 


+0.136 


+0.096 


-0.082 








+0.056 


-0.064 


+0.092 


Pel 


+0.140 


-0.134 


-0.019 


+0.017 


-0.040 


+0.043 


-0.016 


+0.026 


-0.022 


Fell 





+0.051 


+0.084 


-0.067 


-0.023 


+0.021 


+0.034 


-0.034 


+0.063 



Iron and a-abundances in Carina 



Table 5 

Abundance results for the analysed Carina red giants 



LG04a.002181 LG04b.004260 LG04c_006477 



Element 


[X/Fc] 


0"EW 


N 


0"tot 


[X/Fe] 




N 


ftot 


[X/Fe] 




N 


ftot 


Oi 


0.39 


0.40 


2 


0.30 


-0.21 


0.05 


1 


0.08 


0.26 


0.24 


2 


0.19 


Nal 


-0.31 


0.39 


2 


0.29 


-0.86 


0.05 


1 


0.10 


-0.12 


0.38 


3 


0.24 


Mgl 


0.05 


0.03 


2 


0.05 


-0.24 


0.13 


2 


0.11 


0.12 


0.18 


2 


0.13 


Sir 


0.03 


0.10 


4 


0.10 


-0.22 


0.13 


2 


0.17 


-0.13 


0.05 


1 


0.14 


Cal 


0.18 


0.04 


18 


0.07 


-0.02 


0.07 


15 


0.06 


0.06 


0.08 


16 


0.06 


Til 


-0.12 


0.04 


32 


0.12 


-0.17 


0.06 


16 


0.09 


0.13 


0.07 


18 


0.09 


Ti II 


-0.05 


0.11 


6 


0.11 


0.00 


0.07 


5 


0.02 


0.05 


0.05 


5 


0.10 


Fei/H 


-1.38 


0.22 


97 


0.14 


-1.49 


0.16 


96 


0.18 


-1.40 


0.16 


97 


0.18 


Feil/H 


-1.38 


0.18 


7 


0.27 


-1.48 


0.13 


4 


0.20 


-1.37 


0.12 


6 


0.19 




LG04a.000377 






LG04a.001556 






LG04a.OO1826 




Element 


[X/Fe] 




N 


Ctot 


[X/Fe] 


cew 


N 


Ctot 


[X/Fe] 


cew 


N 


Ctot 


Ol 


0.65 


0.50 


■1 


0.36 


0.00 


0.20 


2 


0.16 










Nal 


-0.30 


0.40 


2 


0.29 


-0.55 


0.21 


2 


0.17 


-0.48 


0.06 


1 


0.11 


Mgl 


0.08 


0.19 


2 


0.14 


-0.28 


0.10 


2 


0.08 


-0.90 


0.06 


1 


0.09 


Si I 


-0.07 


0.06 


1 


0.14 


0.14 


0.17 


2 


0.15 


0.06 


0.28 


2 


0.25 


Cai 


0.21 


0.04 


17 


0.04 


0.08 


0.07 


16 


0.04 


-0.16 


0.05 


15 


0.07 


Til 


0.07 


0.06 


26 


0.08 


-0.22 


0.05 


28 


0.20 


-0.83 


0.08 


16 


0.12 


Ti II 


-0.08 


0.07 


5 


0.23 


-0.17 


0.13 


6 


0.11 


-0.63 


0.11 


5 


0.08 


Fei/H 


-1.50 


0.16 


91 


0.18 


-1.27 


0.14 


98 


0.15 


-1.49 


0.16 


95 


0.18 


Feii/H 


-1.49 


0.12 


5 


0.20 


-1.29 


0.14 


7 


0.23 


-1.43 


0.14 


7 


0.22 




LG04C.000777 






LG04d.006628 






LG04C.000951 




Element 


[X/Fc] 




N 


Ttot 


[X/Fe] 


cew 


N 


Ctot 


[X/Fe] 


fEW 


N 


ftot 


Oi 


0.30 


0.08 


1 


0.10 


0.10 


0.06 


1 


0.08 










Nal 


-0.32 


0.13 


2 


0.12 


















Mgl 


-0.03 


0.20 


2 


0.08 


-0.06 


0.13 


1 


0.16 


0.31 


0.12 


1 


0.15 


Si I 


0.01 


0.24 


3 


0.10 


0.11 


0.04 


1 


0.06 










Cai 


0.33 


0.06 


17 


0.09 


0.13 


0.10 


11 


0.13 


0.71 


0.05 


9 


0.07 


Til 


-0.09 


0.06 


29 


0.08 


-0.25 


0.08 


11 


0.14 


0.10 


0.14 


8 


0.09 


Ti II 


0.17 


0.15 


5 


0.15 


0.18 


0.17 


4 


0.14 


0.04 


0.12 


3 


0.10 


Fei/H 


-1.36 


0.12 


90 


0.13 


-1.79 


0.15 


86 


0.18 


-2.72 


0.16 


67 


0.22 


Feii/H 


-1.35 


0.19 


4 


0.29 


-1.79 


0.12 


7 


0.14 


-2.72 


0.12 


4 


0.13 



LG04C-00626 

Element [X/Fe] <tew N CTtot 



Oi 


0.58 


0.05 


1 


0.07 


Nai 


0.02 


0.05 


1 


0.10 


Mgl 


0.41 


0.09 


2 


0.08 


Si I 


0.92 


0.13 


2 


0.18 


Cai 


0.09 


0.11 


9 


0.07 


Til 


0.23 


0.05 


9 


0.10 


Ti II 


0.40 


0.08 


3 


0.03 


Fel/H 


-2.50 


0.15 


79 


0.19 


Feil/H 


-2.54 


0.11 


4 


0.16 



